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We review some recent proposals for relativistic models of dark matter in the context of 
bimetric gravity. The aim is to solve the problems of cold dark matter (CDM) at galactic 
scales, and to reproduce the phenomenology of the modified Newtonian dynamics (MOND), 
while still being in agreement with the standard cosmological model A-CDM at large scales. 

In this context a promising alternative is dipolar dark matter (DDM) in which two different 
species of dark matter particles are separately coupled to the two metrics of bigravity and are 
linked together by an internal vector field. The phenomenology of MOND then results from a 
mechanism of gravitational polarization. Probably the best formulation of the model is within 
the framework of recently developed massive bigravity theories. Then the gravitational sector 
of the model is safe by construction, but a ghostly degree of freedom in the deconpling limit is 
still present in the dark matter sector. Future work should analyze the cosmological solutions 
of the model and check the post-Newtonian parameters in the solar system. 


1 Introduction and motivation 

1.1 Problem of dark matter at galactic scales 

The standard model of cosmology A-CDM is widely held to be an excellent description of reality 
at large cosmological scales. Impressive observational successes of this model include the fit of 
the anisotropies of the cosmic microwave background (CMB), the baryon acoustic oscillations, 
the formation of large scale structures and the accelerated expansion of the Universe. However 
some fundamental issues remain: (i) The measured value of the cosmological constant A looks 
unnatural from a quantum field perspective; (ii) The weakly interacting particles envisaged as 
candidates for the cold dark matter (CDM) are still undetected in the laboratory; (hi) The model 
A-CDM falls short in explaining the observed regularities in the properties of dark matter halos 
around galaxies. 

Regarding point (hi) we have in mind the baryonic Tully-Fisher relation between the observed 
luminous mass and the asymptotic rotation velocity of spiral galaxies the analogous Faber- 
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Jackson relation for elliptical galaxies , and, very important, the tight correlation between the 
mass discrepancy (i.e. the presence of dark matter) and the involved acceleration scale In the 
prevailing view‘d, these issues should be resolved once we understand the complicated baryonic 
processes (e.g. supernova winds and outflows from a central supermassive black hole) that affect 
galaxy formation and evolution. Note that with such foreseen explanation, conclusive tests of 


the A-CDM model become difficult since any new failure can be attributed to some further 
unknown aspect of baryonic physics. More importantly, this explanation is challenged by the 
fact that galactic data are — mysteriously enough — in excellent agreement with the MOND 
(Modified Newtonian Dynamics) empirical non relativistic formula 

A relativistic MOND theory is required to address issues concerning cosmology and grav¬ 
itational leasing. Many relativistic extensions for MOND have been proposed, including a 
tensor-vector-scalar (TeVeS) theory a bimetric theory non-canonical Einstein-iEther the- 

ories ^ Galileon theory a Khronon theory a modified dark matter theory 

Most theories have difficulties at reproducing the cosmological observations, notably the full 
spectrum of CMB anisotropies. An exception is the modified dark matter approach 
which agrees with the model A-CDM at first order cosmological perturbations. This approach, 
also called dipolar dark matter (DDM), is motivated by the dielectric analogy of MOND — 
that MOND represents the gravitational analogue of the Gauss law of electrostatics in dielectric 
non-linear materials ^ natural formulation of the model is based on a bimetric extension 

of general relativity (GR) where two species of dark matter particles are respectively coupled 
to the two metrics, and are linked by an internal vector field generated by the mass of these 
particles. In this model the phenomenology of MOND emerges naturally and elegantly from a 
mechanism of gravitational polarization. 

1.2 Bimetric massive gravity theories 

Bimetric theories have been extensively investigated in the quest of a consistent massive gravity 
theory. From a theoretical point of view, the existence of a graviton mass is a very important 
fundamental question. At the linear level, there is a unique mass term which ensures the absence 
of ghosts in the theory, namely the Fierz-Pauli action . Even though theoretically consistent, 
this linear theory suffers from the van Dam-Veltman-Zakharov (vDVZ) discontinuity which 
reflects the fact that one does not recover GR in the limit of vanishing graviton mass. In 
other words, this discontinuity can be traced to the coupling of the additional longitudinal 
graviton to the matter field in the mass going to zero limit. Vainshtein very soon realized 
that the nonlinearities of the theory become actually stronger in the vanishing mass limit , 
pointing that these nonlinearities might cure the vDVZ discontinuity, which required the non¬ 
linear completion of the Fierz-Pauli theory. Unfortunately, this task seemed to face the inevitable 
problem of reintroducing the ghost-like instability^^. 

Recent joint effort to construct a consistent ghost-free non-linear theory for massive gravity 
gave fruitful results which have initiated a revival of interests. The theory can be 

further generalized to bigravity or multigravity , by the inclusion of the corresponding 
additional kinetic terms. Another fundamental question in the context of massive gravity is 
the consistent coupling to the matter fields. If one couples the matter fields in a naive way 
additively to both metrics simultaneously, this immediately reintroduces the Boulware-Deser 
(BD) ghost Moreover, quantum corrections at one loop dictate that the only consistent 

way of coupling the matter sector to the two metrics is either through a minimal coupling to 
just one metric (which preserves the technical naturalness of the theoryor through a 
composite effective metric built out of both metrics in a very specific way'^^’'^^’'^'. An important 
consequence of the coupling through the effective metric is a possible way out of the no-go 
result for the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric together with the 
propagation of the five physical degrees of freedom of the graviton sector without introducing 
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ghost and gradient instabilities'^®. 

Massive gravity is replete of phenomenology. Especially, its potential application in cosmol¬ 
ogy received much attention. Even if the decoupling limit of the theory admits self-accelerating 
solutions , the full theory with Minkowski reference metric suffers from the no-go result for 
flat FLRW solutionThe cosmology of the bigravity theory has more freedom and features 
due to the dynamics of the reference metric Assuming that the matter fields couple 


minimally to one metric and further assuming that the mass of the graviton is of the same order 
as the Hubble parameter today, the theory admits several interesting branches of solutions. 

1.3 Content and summary 

In Sec. 2 of this paper we shall review the initial bimetric model ^ and the broad range of 
phenomenology it is able to predict. On the bad side of this model, we shall also discuss, in 
Sec. 2.3, with the help notably of the minisuperspace of the model, the likely presence of ghost 
instabilities. This will motivate the redehnition, in Sec. 3, of the gravitational sector of the 
model ^ in a way to make it consistent with the beautiful framework of massive bigravity 
theoriesGladly, this move will substantially simplify the model. On the other hand, the 
dark matter sector will essentially remain the same as in ^. The mechanism of gravitational 
polarization, and recovery of the MOND equation, checked in Sec. 3.2, thus appear as a natural 
consequence of massive bigravity theory for this type of matter. Concerning ghosts, the new 
model is safe in the gravitational sector (by construction), but a ghostly degree of freedom in 
the decoupling limit is reintroduced in the dark matter sector . A crucial question to address 
in future work is whether the polarization mechanism can be realized in absence of ghosts. 


2 Bimetric theory with two dark matter species 


2.1 Relativistic action 


A relativistic model involving, in addition to the ordinary matter simply described by baryons, 
two species of dark matter particles, was proposed in^*^. A vector field is sourced by the mass 
currents of dark matter and dubbed graviphoton. This vector field is crucial in order to ensure 
the stability of the dipolar medium. The gravitational sector is composed of two dynamical 
Lorentzian metrics and The baryons (representing in fact the full standard model of 
particle physics) are coupled in the usual way to the metric The two species of dark matter 
particles are respectively coupled to the two metrics g^^ and . The gravitational-plus-matter 
action of the model reads (in geometrical units G = c = 1) 
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Here Rg, Rf and TZes are the Ricci scalars of the g and / metrics, and of an effective composite 
metric defined non-perturbatively from g^i, and fgi, by“ 


Gf. = = UpYS , 


( 2 ) 


where the square root matrix is defined by A = g together with its inverse Y = / ^g. 

Notice that Gf^ can be shown to be automatically symmetricThe dimensionless coupling 
constant e measures the strength of the interaction between the two sectors g and /. We inserted 
three different cosmological constants Xg, A/ and Aefr in the g, f and interaction sectors. They 
will be ultimately related to the observed cosmological constant A. 

Baryons and dark matter particles are described by pressureless fluids with conserved scalar 
densities pbar, Pg and pf. In addition Jg and stand for the mass currents of the two types 

“Actually the effective composite metric was defined in " by the different condition 

_ ^P<7 p _ ^pCT p 

yfiv — yefldpi^ Jv<T — Ueff • 

It has been shownthat this condition is equivalent to the requirement (2). 











of dark matter particles, defined by 

, ( 3 ) 

where Jg = pgUg and Jj = Pfu'^ are the conserved dark matter currents associated with the 

respective metrics g and /, thus obeying '^fiJg = 0 and = 0. 

The vector field Ag obeys a non-canonical kinetic term VV{X) where 

with Fgy = dgAu — dyAg. Note that the vector held strength in (4) has been rescaled by the 
MOND acceleration oq — 1-2 10“^*^ m/s^. The function W’(T) is determined phenomenologically, 
but in principle it should be interpreted within some more fundamental theory. In the limit 
T ^ 1, which corresponds to the MOND weak-acceleration regime below oq, we impose 

>V(T) = T - + O (T^) . (5) 
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On the other hand we also impose that when T S> 1, corresponding to the strong-acceleration 
regime much above uq, 

>V(T) = A + -^ + 

where A and B are constants and b > 0. 


2.2 Phenomenological predictions 

We now review the rich phenomenology of this model, which is quite successful at different scales 
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and in different regimes. We refer to ^ for the details. 


• Cosmology. We study the cosmology of the model by expanding the two metrics around 
two homogeneous and isotropic FLRW background metrics. 


dSg = Og [—dr/^ -|- dx* dx-^] , (7a) 

dsj = Oj [—dr/^ -|- dx* dx-^] , (7b) 


where g is the conformal time and Og and ap are the two different scale factors. We show 
that the consistency of the equations in the background, where the two metrics superpose, 
is ensured provided that the background densities obey. 


Pb 
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where a = ag/ap is constant, and p is the common density of the two dark matter fluids 
in the background. The cosmological constants should also be related to the observed 
cosmological constant A, through the relations. 


Xg = A, ~ o? A , Aeff = a. A. 


(9) 


Then, at first order cosmological perturbations around the FLRW metrics, we define the 
g'-sector as being the observable one because it is where the baryons live and light propa¬ 
gates, and all observations take place. We thus study the perturbation equations in this 
sector. We define some effective dark matter quantities that modify the ^-type dark mat¬ 
ter particles taking into account their interaction with the other sector. Using these new 
variables we find that the perturbation equations in the g sector take exactly the same 
form as those of the A-CDM model. Thus the model turns out to be indistinguishable 
from A-CDM up to first order perturbation, and is then fully consistent with the observed 
fluctuations of the CMB. 






• MOND. At galactic scales, in a regime of small accelerations a <C QQ) the potential function 
W takes the form (5). Based on a particular solution for the equation of the vector held 
Afj,, the two dark matter huids can be described as a polarizable dipolar medium. The 
MOND equation is then recovered as a result of a mechanism of gravitational polarization 
which appears as a natural consequence of the model. Moreover the dipolar dark matter 
medium undergoes stable plasma-like oscillations. We shall give some more details of this 
polarization mechanism in Sec. 3.2 at the occasion of the next model based on massive 
bigravity. However note that the next model has essentially the same predictions with 
regards to MOND phenomenology as the model There is though a difference, in that 
in the model we have to assume that the coupling constant e in the action ( 1 ) tends to 
zero, while no particular requirement will be necessary in the next model. 

• Solar System. As we modify GR we have to check the post-Newtonian limit in the Solar 

System, for which the potential function W in the action takes the form ( 6 ). We have 
shown that when we expand the model at first post-Newtonian order, the parametrized 
post-Newtonian (PPN) parameters are exactly the same as in GR. A crucial point for 
this test is a non-linear effect present in our definition of the effective metric ( 2 ) and 
which permits to recover the correct value for the parameter measuring the amount of 
non-linearity, = 1 . Thus the model passes the Solar System tests and is viable. 


2.3 Minisuperspace and ghosts 
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Unfortunately, the previous model proposed in ^ suffers from an harmful ghost instability. The 
source of its origin is multifaceted. First of all, the presence of the square root of the determinant 
of GeS in the action ( 1 ) corresponds to ghostly potential interactions, 


V~GeS = \l■ ( 10 ) 

Already at the linear order around a flat background, posing and = 

{iliiu + these potential interactions do not preserve the Fierz-Pauli tuning, with [• • • ] 

denoting the trace as usual, 

7=5^= 1 + i[(. + f] + i([ft + + f2]) 

+ + + + + +■■■ , ( 11 ) 

yielding a linear ghost instability already present at the scale 
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where A 3 = (Mpim^)^/^ and vr encodes the helicity -0 degree of freedom of the massive graviton 
with mass m. The ghost corresponds to a very light degree of freedom, and hence the theory 
cannot be used as an effective field theory. The other source of ghostly interactions in the 
model is originated in the presence of three kinetic terms. Studying the theory in the mini¬ 
superspace immediately reveals the ghostly interactions of the considered kinetic terms. The 
dangerous sub-Lagrangian due to kinetic interactions is given by 
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As a first diagnostic, we can assume the two metrics to be of the mini-superspace form 


dSg = 5/2i/dx^dx^ = —Ugdt'^ -|- a^dx^ , 


(14a) 














= —n'fdt'^ + ajdx'^ , 


(14b) 


where Ug, Uf are the lapse functions and ag, aj are the scale factors of the respective metrics. 
The effective metric Qes in the mini-superspace becomes 

+ algdx"^ , (15) 


where the effective lapse and effective scale factor correspond to rieg = ^Jngnf and Oeff = ,Ja^aJ. 
We compute the conjugate momenta for the scale factors and get 


Pg — QM-p^agHg [HgUg + HfUf), 

4 UeS 

Pf = - ^M^sag^{Hgng + HfUf) , (16) 

^ npff 


with the conformal Hubble factors Hg = and Hj = respectively. After performing the 
Legendre transformation we obtain the Hamiltonian in the mini-superspace: 


eff ^ (PgQg - Pfaf?ng 
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As one can see immediately, the Hamiltonian is highly non-linear in the lapses Ug and Uf, 
signalling the presence of the BD ghost degree of freedom already in the mini-superspace. Sum- 

on 

marising, the model proposed in^^ cannot represent a consistent theory for dipolar dark matter 
in bigravity due to the ghostly contribution in form of the cosmological constant for and 
the kinetic term V~GesT^eS- 


3 Model based on massive bigravity theory 


3.1 Covariant theory 
on 

The previous model ^ is plagued by harmful ghosts, but it remains that the phenomenology, 
especially at galactic scales (i.e. MOND), calls for a more fundamental theory. We now look 
for a consistent coupling of the dark matter particles within the framework of massive bigravity 
theory and the restrictions made in 32,36,37 concerning matter helds. We thus propose a new 

on 

model, whose dark matter sector is essentially the same as in the previous model ^, but whose 
gravitational sector is now based on ghost-free massive bigravity theory. As bigravity theory 
represents essentially a unique consistent deformation of GR, we think that the new model could 
represent an important step toward a more fundamental theory of dark matter a la MOND in 
galactic scales. 

The model is based on the bigravity-plus-matter action 


S = 



(18) 


where Mg and Mf are two coupling constants, and m is the mass of the graviton. The ghost-free 
potential interactions between the two metrics g and / take the particular form of the square 

on AQ 

root of the determinant of the effective metric 


9‘^!^ = + 2a(3 -h (3“^ , (19) 

where a and (3 are arbitrary constants, which can always be restricted in the model (18) to 
satisfy a + f3 = 1, and denotes the effective metric of the previous model as defined by (2) 












in terms oi X = sj g ^ f and Y = The square root of the determinant of this effective 

metric admits a closed-form expression, 


= diet[al +fix) = ^/^ <\ei[l31 + aY) , (20) 

and can also be written with the help of the usual elementary symmetric polynomials CniX) 
and en{Y) of the square root matrices X or T as 

4 4 

= «"/5"""en(T). (21) 

n=0 n=0 


We see that (20)-(21) corresponds to the right form of the acceptable potential interactions 
between the metrics g and /. To recover the usual Newtonian limit for the motion of baryons 
with respect to the ordinary metric g we must impose 


Svr 


( 22 ) 


Finally the vector field Afj, is now coupled to the metric 5 ®® rather than to thus 


X = - 



2aQ 


(23) 


However the functional form of the non-canonical kinetic term W in the MOND regime is the 
same as in the previous model, see (5), and we could also impose ( 6 ). Finally the mass currents 
to which is coupled the vector field in (18) are dehned by 


V-5efrj^ = , y/-9esjf = > 


(24) 


where, as before, Jg = pgUg and Jj = PfU^- However, notice the presence of the factor a/fi in 
the interaction term between Ag and the current in (18). This factor can be interpreted as 
a ratio between the gravitational charge and the inertial mass of the / particles when measured 
with respect to the g metric. 


3.2 Gravitational polarization & MOND 

We now discuss the main point of this model, which is to allow a mechanism of gravitational 
polarization that permits to recover in a straightforward way the phenomenology of dark matter 
at galactic scales (MOND). In this respect the predictions of the new model are essentially the 
same as those of the previous model ^. 

For this purpose we are mainly interested to that particular combination of the two metrics 

zl.7 

g and / which is massless '. Working in the non-relativistic limit c 00 , we find that the 
massless combination reduces to an ordinary Poisson equation for the Newtonian potentials Ug 
and Uf associated with the two metrics, namely 

( 2M2 2M2 \ 1 1 

where /O^ar’ P*g P*f denote the ordinary Newtonian densities of the matter fluids. However, 
with massive (bi-)gravity the two sectors associated with the metrics g and / do not evolve 
independently but are linked together by a constraint equation which comes from the Bianchi 
identities. We find that this constraint reduces in the non-relativistic limit to 


V{aUg + (iUf) = 0 . 


( 26 ) 











Combining (25)-(26) and using (22) we readily obtain the following Poisson equation for the 
Newtonian potential in the ordinary sector, 


AUg = -A7:[pl,, + pl--^p}). 


(27) 


Similarly we find that the equation governing the vector field, obtained by varying the action (18) 
with respect to A^, reduces in the non-relativistic limit to a single Coulomb type equation, 
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(28) 


where Wx is the derivative of W with respect to its argument X. Finally we have also at 
our disposal the equations of motion of the baryons and the dark matter particles. We look for 
explicit solutions of the equations of motion of dark matter in the form of plasma-like oscillations 
around some equilibrium configuration. Thanks to the constraint equation (26) we find that an 
equilibrium is possible, i.e. the dark matter particles are at rest, when the Coulomb force 
annihilates the gravitational force, namely 


Vt/g -|- V(/> = 0 . 


(29) 


At equilibrium we shall grasp a mechanism of gravitational polarization, i.e. we can define a 
polarization field which will be aligned with the gravitational field. Away from equilibrium 
we find that the polarization field undergoes stable plasma like oscillations (see^^ for details). 
Finally, combining the three equations (27), (28) and (29) we easily deduce that the equation 
for the ordinary potential Ug takes exactly the Bekenstein-Milgrom form'^^. 
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(30) 


with MOND interpolating function p = 1 —Wx- Furthermore, thanks to the postulated form (5) 
of the function yV[X) we recover the correct deep MOND regime when X —)■ 0. We refer 
to 19-20,45 £qj. j] 2 ore details about this way of recovering the MOND phenomenology, which is of 
course reminiscent of the dielectric analogy of MOND. 

Let us emphasize that gravitational polarization &: MOND appear as natural consequences 
of this model, and are made possible by the constraint equation (26) linking together the two 
metrics of massive bigravity theory. Furthermore this requires that the gravitational force can 
be annihilated by some internal non-gravitational force, here chosen to be a vector field. This 
implies the existence of a coupling between the two species of dark matter particles living in 
the g and / sectors. Unfortunately, as has been shown in'^^, the latter coupling is problematic 
as it yields a ghostly degree of freedom in the dark matter sector. Further work is needed to 
determine the exact mass of the ghost and see whether the required polarization mechanism and 
the ghost absence are compatible. On the phenomenology side, the cosmological implications 
of the model and the post-Newtonian parameters in the Solar System should be investigated in 
great details, as it has been done for the previous model 
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